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Without involving bounce events, a Poincare section associated with the axes is found to give 
a map on the annulus for the diamagnetic Kepler problem. Symbolic dynamics is then established 
based on the lift of the annulus map. The correspondence between the coding derived from this axis 
Poincare section is compared with the coding based on bounces. Symmetry is used to reduce the 
symbolic dynamics. By means of symbolic dynamics the admissibility of periodic orbits is analyzed, 
and the symmetry of orbits discussed. 
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I. INTRODUCTION 



The diamagnetic Kepler problem (DKP), or a hydrogen atom moving in a uniform magnetic field, is a simple 
physical system for studying the correspondence relation between classical chaos and quantum behavior ffl. In the 
semiclassical quantization of a classically chaotic system, quantum properties are related to periodic orbits or closed 
orbits An effective method to locate orbits, and then encode and classify them is essential to the study. For the 

n-disk scattering system || and stadium billiard symbolic dynamics has been established by taking a Poincare 
section on the boundary of reflection. The symbolic encoding of the orbits and symbolic dynamics developed for the 
four-disk system have been extended to the smooth Hamiltonian system for the DKP [|s|-[l0||. However, bounces are 
^« . not always well defined for a soft potential. For the hyperbolic potential system, symbolic dynamics is established 
by introducing an appropriate Poincare section without referring to bounce events Here we use this method to 
OS , construct symbolic dynamics for the DKP. The paper is organized as follows. In Sec. [IlL we describe the method to 
construct a Poincare section with the Birkhoff coordinates. In Sec. |Hl| we discuss the symbolic encoding of orbits 
and ordering rules of symbolic sequences. We examine the correspondence between this coding and the coding based 
£>y on bounces for the touching four-disk billiard in Sec. IV. In Sec. |v|, we classify orbits according to their codes and 



symmetries. Finally, in Sec. VI, we give a summary and conclusion. 



o 

i II. THE POINCARE SECTION 

o 

The classical dynamics of a hydrogen atom in a uniform magnetic field B along the z-axis is described by the 
1 Hamiltonian: 

H = ^--+ujl z + ]-m e oj 2 {x 2 +y 2 ) (1) 
2m e r 2 

where l z is the z component of the angular momentum, and u> = -^- c the cyclotron frequency. In cylindrical coordinates 
(p,z,4>) and atomic units, at the angular momentum l z = 0, the Hamiltonian can be rewritten as 

H -Y + T~Y-^T7i> (2) 

which is an even function of the coordinates and momenta. Its symmetry group consists of the identity e and the 
reflection a p . Introducing the semi-parabolic coordinates 



/i 2 = vV + z 2 



v 2 = \/p 2 + ^ 2 + z, (3) 

(4) 
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we may convert the Hamiltonian to 



where e = E(2u>) 2 / 3 is the scaled energy and /i is the fixed "pseudo energy". In Fig. 1, we show an orbit at e = 0. 
Irregular orbits at e = displayed in the configuration space(z, p) Q or in the Poincare section [v, p v ) show the 
chaotic behavior of the dynamics. 

The transformed Hamiltonian h has a higher symmetry than H. Its symmetry group consists of the identity e, two 
reflections er M , a v across the /i, v axes, two diagonal reflections <7i3, CT24, and three rotations C4, C2 and Cf by 7r/2, 



7t and 37r/2 around the center, respectively 12 



Our Poincare section is chosen as follows. Imagine that the fj, and v axes are both of a finite width and length. A 
counterclockwise contour is taken along the perimeter of the area forming by the two crossing imaginary rectangles. 
The Poincare section is then obtained by recording the position and the tangent component of the momentum along 
the contour, i.e., the Birkhoff canonical coordinates |l3f at intersecting points with the contour where an orbit enter 
the inside of the contour. The length of the contour is infinite. It is more convenient to transform the contour to one 
with a finite length. For example, in the first quadrant, the transformations s = —fj,/ (1 + /1) along the positive /x axis 
and s = + v) along the positive v axis convert the segment of the original contour in the first quadrant into an 
interval of length 2 parametrized with s G [—1, 1). For our purpose here, the variable corresponding to the momentum 



may be take as v = —p^/p at the positive fi axis, and v = p v /p at the positive v axis, where p = yP 2 +pi- I n this 

way we may parametrize the whole contour with s G [—1,7), and define corresponding v. The rotational symmetry 
under C4, C2 and C| in the original configurational space becomes the translational symmetry of shifting s by a 
multiple of 2 in the s-v plane. The dynamics on the Poincare surface is then represented by a map on the annulus 
s G [-1, 7) and v G [-1, 1]. 



III. SYMBOLIC DYNAMICS BASED ON LIFTING 



For the map on an annulus describing the dynamics on the Poincare section, it is useful to consider its lift. In the 
lifted space on which the lift map is defined the area s G [—1,7) and v G [—1,1] consists of a fundamental domain 
(FD), which is uniformly divided into 8 strips according to s. The image of the FD in the lifted space is partly shown 
in Fig. |^(a), where the 8 strips of the FD are marked with the numbers to 7, and four strips right of the FD marked 
with +0 to +3. In the figure we draw only the images of the strips 1 and 2. Two zones marked with 1 and 1' of the 
strip 1 are mapped to the zone 1 of strip +3 and the zone 1' of strip +2, respectively. The notation representing the 
image of strip 2 is analogous. The rotational symmetry of the Hamiltonian corresponds to the translational symmetry 
in the lifted space. For example, the image of strip is in the strips 6 and 7, and can be obtained from the image of 
strip 2 (in the strips +0 and +1) by shifting to the left by 2 strips. Similarly, we show part of the preimage of the FD 
in Fig. |2|(b). Since on the annulus forward and backward maps rotate in opposite directions, strips are on the right 
of their preimages, in contrast with images. 

For the full FD, the rotation number of orbits is between and 2. Taking the rotational symmetry into account, we 
may consider only the reduced domain (RD) consisting of strips 1 and 2. When the RD is regarded as the annulus, 
from Fig. 1 the integral part of the rotation numbers is 3 and 4 for zone 2, which is the joint zone of 1' and 2', but 4 
and 5 for zone 1. As shown in Fig. ||(a), the RD may be correspondingly divided into three regions according to the 
rotation number. The two lines of demarcation are marked by Bq» and i?2»- In the figure we have also drawn some 
stable and unstable manifolds. (The stable manifolds have a segment roughly parallel to lines Bqu and B2*, and the 
unstable manifolds are a mirror image of the stable manifolds with respect to the line v = 0.) In the RD, tangencies 
between stable and unstable manifolds can be easily seen. Among them, the prominent ones are in the top-left and 
bottom-right regions of the three. At a given point the stable and unstable directions, i.e. tangent directions of 
manifolds, can be determined with the procedure suggested by Greene [ [Ti) , and then tangent points found. Two 
lines connecting such tangent points, which are marked by Co* and C 2 *, give a further partition of the RD. The final 
partition of the RD into five regions is shown in Fig. ||(a), where these regions are marked by L », R », Ri», R 2 » 
and L2*. In the figure the two boundaries of the RD are marked by Dq» and D 2 *. This partition is the partition 
according to the preimage. Its preimage gives the partition according to the present, as shown in Fig. ||(b), where the 
five regions are marked with »Lq, »Rq, *Ri, »i?2, and *Li- By means of this partition, or equivalently the partition 
of Fig. [|(a) according to the preimage, we may code an orbit with a doubly infinite sequence 

• • ■ S_l • SqSi 
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where • indicates the present. We may refer to this coding as the axis coding, and the coding based on bounces in 
the literature the bounce coding. 

The ordering is essential to establishing symbolic dynamics. The natural order is well defined in the lifted space. 
The 5-piece preimage of the RD arranged from left to right in the lifted space consist of Lqu, Rq*, Ri*, R 2 * and L 2 u, 
as can be seen from Fig. ||(b). Thus, according to the natural order we have 

L » < i? » < Ri» < i?2» < L 2 • . (6) 

From the preimage of the RD shown in Fig. |^(b) we see an opposite arrangement of the preimage pieces in the lifted 
space. A more precise description of ordering is given by the ordering of unstable manifolds on a transversal stable 
foliation. It is numerically verified that under the backward map the ordering is reversed in the regions Rq», R\» and 
i?2», but preserved in Lqu and L 2 ». We may define the parity of a finite string by the oddness of the total number 
of the letters Rq, Ri and R- 2 contained in the string. Any odd leading string in backward sequences will then reverse 
the ordering (||) of backward sequences. 

Similarly, from the partition shown in Fig. [|(b) according to the present, we have the ordering 

• L < uR < < »R 2 < »L 2 . (7) 

An odd leading string also reverses this ordering for forward sequences. 

Based on the ordering rules,metric representations for both forward and backward sequences may be introduced to 
construct the symbolic plane [|15| . Every forward sequence »sqSi ■ ■ ■ s n ■ ■ ■ may correspond to a number a, represented 
in base 5, between and 1 as follows. We define the correspondence of symbol Si to number /i^ 6 {0,1,2,3,4} as 
Lq — > 0, Rq — > 1, R\ — > 2, R 2 — > 3 and L 2 — > 4 if the leading string sqSi ■ ■ ■ Sj_i is even, and as {Lq, Rq, Ri, R2, L 2 } — > 
{4, 3, 2, 1, 0} otherwise. Finally, we define 



i=0 



(8) 



In this way forward sequences are ordered according to their a-values. Similarly, we assign to a backward sequence 
• • • s_ m • • • s_ 2 s-i» the number j3 defined by 

oc 

= £"i5 _i , ( 9 ) 
3=1 

where i/j 6 {0, 1, 2, 3, 4} is determined by s_j and the oddness of the leading string S-j+i • • • S-2 s -i- In the symbolic 
plane every orbit point corresponds to a point (a, (3) in the unit square, where a and f3 are associated with the 
forward and backward sequences of the orbit point, respectively. We show the symbolic plane for e = in Fig. |](a), 
where approximately 68000 points of several real orbits in the RD are drawn. The corresponding pruning front for 
the partition lines «Cq, »C 2 and the boundary lines »D and »D 2 is shown in Fig. |](b). 

So far we have not considered the reflectional symmetry. In the lifted space this symmetry corresponds to the 
invariance of the RD under a 7r- rotation around the center (or the reflection with respect to the center). From Fig. 
^(a) it is seen that the image of zone 1 of strip 1 is in strip +3, hence still in 1 after wrapping. On the contrary, the 
image V of strip 1 is in strip 2. However, a 7r-rotation can put the image back into strip 1. So, using the reflectional 
symmetry, we may focus only on strip 1, which may be regarded as the minimal domain (MD). In this way the 
5-letter symbolic dynamics is reduced to a 3-letter one, and the 7r-rotation changes the parity of Ri . More specifically 
speaking, we have the following ordering for the minimal domain 

• L < »Ro < •Ri, L » < Rq* < Ri», (10) 

and only Rq in a leading string reverses the ordering. The symbolic plane of the MD is shown in Fig. ^ where 
34,000 points of several real orbits are drawn together with the pruning fronts. The primary pruning front of »Dq 
is roughly diagonal, while that of »Co is almost vertical. The latter encloses a forbidden zone near the top in the 
symbolic plane, which is responsible for many interesting bifurcations. In the 5-letter code for the RD, the second 
and fourth quadrants of the symbolic plane are forbidden. The symbolic plane for the 3-letter code of the MD looks 
more compact. 

Let us discuss the pruning front in some detail. A point Q • DqP on the partition line »Dq discards a rectangle 
in the symbolic plane as a forbidden zone, the top-right corner of which is the point Q • DqP, as shown in Fig. |^. 
When representing the point in the symbolic plane, we have replaced Dq by its right limit, i.e., a point with s = 0+. 
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Its allowed zone is the rectangle, the lower-left corner of which is the same point. Similarly, a point Q • CqP discards 
a rectangle as its forbidden zone whose two lower corners are Q • LqP and Q • RqP. The point also determines 
two allowed zones. The above mentioned two corners form the upper-right corner of the left allowed zone and the 
upper-left corner of the right one. We have numerically found the following 5 points on the partition lines: 

Pi '. ■ ■ ■ R\L^R\LqR,\LqR\R^R\ • CoR\RoR\R^R\R^R\ • • • , 

P% \ ■ ■ ■ RiLqRiRqRiLqRqLqR^ • CqRiRqRiRqRiLqRiRqRi • ■ • , 

P 3 : ■ ■ ■ RIL R,IrILqRiRoRi • CqR\RqR\ 2 L • ■ • , 

P4 : • ■ ■ RiLqRiLqRiLqR-^RqRiLqRiLqRq • LqRiLqRiLqRiRqR-^LqRiLqRi • • • , 

-P5 : • • • R\LqR\LqR\RqR\LqR\LqR\Lq • RqPqRiPqRiRqR^LqRiIjqRi • • • , 

where the last two points are on »Dq. The first tree points correspond to lower corners of the top three rectangles with 
dotted sides in Fig. |L while the other two to the upper-right corners of the two bottom rectangles. It can be verified 
that the point (RqRiRqRiLqR\)°° • (RqR^RqRiLqR^) 00 is in the forbidden zone of P2, so the periodic sequence 
(RoRiRqRiLqRI) 00 is nonadmissible. One can verify the admissibility of the sequence (R\LqRiRqR\R^RiRqRi)°° . 
In fact, its 0, 1, 2, 4, 5, 7, 9 th shifts are in the allowed zone of Pi, while its 3, 8 th shifts in that of P3. At the same 
time, its 0, 1, 3, 5, 6, 8 th shifts are in the allowed zone of P4, while its 2, 4 7, 9 th shifts in that of P3. In this way 
we may examine the admissibility of periodic sequences from the known coding of partition lines. The above orbit 
(RiL RiR R\RqRiR q R 1 )°° , written in the 5-letter code, is {RiR2RiRoR\RqRiL2Ri)°° ■ Symmetry of sequences by 
means of the 5-letter code will be discussed in Sec. V. 



IV. CORRESPONDENCE BETWEEN THE AXIS CODES AND THE BOUNCE CODES 

In the literature the symbolic description based on bounces is used. Bounces are well defined for billiards. We 
may explore the relation between the axis codes and the bounce codes by examining a billiard. For simplicity, let us 
consider the touching four-disk billiard. As shown in Fig. 0(a), the boundary of the billiard is formed by the four 
disks which are labeled with a, b, c, and d, respectively. Taking the boundary as a Poincare section, we have the 
corresponding Birkhoff coordinates similar to those of the stadium billiard, which we denote by (r, u) to distinguish 
from the above (s,v), and hence the associated map on the annulus JtJ. We normalize the boundary length as unity, 
and measure the length from the touching point of disks a and d. As shown in Fig. 0(b), the fundamental domain 
of the annulus consists of four strips corresponding to the four disks. The image of the strip of disk a, which has 
three pieces in three other strips, determines the partition of the reduced domain consisting of a single strip. The 
three regions of the RD are labeled with »0, »1, and »2 in the figure. By the argument based on lifting, we have the 
ordering 

• 0<»1<«2. (11) 

Since every bounce on a disk reverses the ordering of orbits, each symbol has an odd parity. That is, for forward 
symbolic sequences, a string of an odd length is odd, while that of an even length is even. We have similar ordering 
rules for backward symbolic sequences. This is our version of the symbolic dynamics described in Ref. 

When the strip of r € [0,0.25) is divided at r = 0.125, and an eighth of the FD is considered, the image of the 
line r = 0.125 has three pieces respectively in strips b, c and d. The partition of a half strip, which is determined by 
putting all the pieces of its image together back into the half strip, will involve symbols more than three. 

To explore the correspondence between the axis codes and the bounce codes, we examine how an orbit starting 
from a point on an axis hits different disks. This gives the relation between the phase spaces (s,v) associated with 
the axis contour and (r, u) associated with the billiard boundary. Consider the orbit AC grazing against disk d in Fig. 
^(a). It corresponds to a point of s 6 ( — 1, 0), or equivalently a point of s £ (1, 2) when the RD is considered. Starting 
from point A on the x-axis, orbits right of this grazing orbit hit disk d before intersecting an axis again, while orbits 
on its left hit the y-axis directly before having a bounce on disk c. The orbit ABO, which passes through the origin, 
separates orbits crossing the y-axis from those crossing the x-axis. The orbits of the type AC correspond to line »Ci 
in the RD of the s-v phase space since the line demarcates the order preserving region without a bounce from the 
order reversing region with a bounce. Similarly, the orbits of type ABO correspond to line »B2- Thus, bisecting the 
region into *R + and with line s = 1, we may re-label the right half of the RD as in Fig. ||(b). Comparing 
this figure with Fig. 3(b), we see that »R+ and »i?2 belong to »0, which associates with orbits leaving from disk a 
and bouncing on disk d. The region »L2 has no direct correspondence with any bounce. Orbits from the region, after 
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missing a bounce, will hit disk c before intersecting again an axis. We may re-label the region as (»1). Consequently, 
we have the correspondence: 

»R+, *R 2 -v «Q; »L 2 X^*1, 

where X stands for any symbol. 

The other half of the RD in s-v space may be analyzed in a similar way. The re-labeling of regions is also shown 
in Fig. |8|(b), which indicates the correspondence: 

•R-,»Ro -» «2; »L Q X -> *1. 

Putting Fig. ||(a) on the top of 3(b), we see that the allowed symbolic pair are 

L2R+, L2R2', R2R+, R2R21 R2L2', R—R+,R—R2,R—Li2; R+R-, R+Rq, R+Lq; RqR-, RqRq, RqLq; LqR_,LqRq. 
Two rules can then be extracted: 

1. L2 and Lq are always followed by a symbol of odd parity. 

2. R+ follows either L2, or R2, or R-, while i?_ follows either Lq, or Ro, or 

The first rule explains that LqX or L2X has the same parity as the bounce code 1. The second rule can be used to 
refine R\ into R + and and then further to convert them into bounce codes. For example, an orbit shown in Fig. 
9(a) has the axis codes (R^RiRqRi) 00 ■ The codes are first refined as {R\R + RqR-)°° , and then converted to (0 3 2 3 )°°. 
Another orbit shown in Fig. 9(b) has the codes (R^RILqRqRi) 00 . This sequence is converted as follows: 

(R^RiLqRqRi) 00 — * {R^R+R—R+LqRqR—) 00 — ► (0 2012) 00 . 

The ordering (0) can be refined as 

•Lq < *Rq < *R- < *R+ < *R 2 < *L 2 , 

which, under a cyclic transformation, becomes 

•R + < *R 2 < *L 2 < *L < »R < mR_. 

This ordering is then converted to 

•0 < »1< »2, 

which is just the ordering of the bounce codes. A detailed analysis of the relation between axis codes and bounce 
codes will be presented elsewhere. 

V. SYMBOLIC SEQUENCES AND SYMMETRY OF ORBITS 

The configurational space of the original dynamics is the half p-z plane, which converts to the first quadrant of the 
p-v plane. The rotational symmetry in the p-v plane is the invcirieincG under the transformation from v } p v ) to 
(—v,fi, —pv,Pn), (— (-1, —v, —p^, —p v ) and (y, —p,p y , —p^) (denoted by p, ir and p). This symmetry has nothing to do 
with the original dynamics. Confining the dynamics to the RD, we remove the symmetry. 

As for the reflectional symmetry, we may focus only on : (p, — v, p^, —p v ) since we can composite a v = n o a^, 
CT13 = 0^0 p and (T24 = &[i p- We shall denote simply by a. While p, n and p have no effect in the s-v space, 
a leads to (s,v) — > (2 — s, —v), which, written in symbols, is Lq <-> L2, Rq <-» R2 and R\ <-» Ry. To discuss the 
reflectional symmetry, the 5-letter code must be used. Thus, if an orbit is symmetric under reflection, its symbolic 
sequence is invariant up to a shift under the symbolic transformation of a. For example, sequence (R\R\R\R^)°° of 
period 6 is symmetric under a, while (R^Rq) 00 is not. This two orbits are shown in Fig. 9(a). 

The Hamiltonian dynamics has the time-reversal symmetry: (p, v, Pi_ t ,p^; r) — > (//, v,P)j,,Pu'i — t) = 
(p, v, —p^, —p v \ t). This transformation, denoted by T, leads to (s, v) — > (2 — s, v) in the s-v space. Written in symbols, 
this means »Lq <-» L2*, *Rq «-> R2*i *R\ Ri*, *R2 <-> Rq* and »L2 <-> Lq». Thus, the effect of T on a sequence 
is to reverse the sequence and apply a to it. For example, transformation T converts sequence (RiLqRqRiRI) 00 to 
(RqRiR2L2R\)°° ■ We may combine T and a to define the transformation T a of reversal by T a = T 00 = o oT . In the 
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3-letter symbolic dynamics we identify L 2 with Lo, and R 2 with Rq, so T just reverses a sequence. Transformation T a 
converts the above just mentioned orbit (Rf LqRqRxRI) 00 to (R%RiRqL RI)°° . The two orbits (Rf LqRoRxRI) 00 and 
{RIRxRqLqRD 00 are shown in Fig. 9(b). (In the configuration space orbits (i?fi -Ro-Ri-Ri)°° and (R 2 l R 1 R 2 L 2 Rf) aD 
are not distinguishable.) 

The transformations a, T and T a may be used to generate admissible sequences from a known one. When performing 
a transformation on a sequence results in the same sequence up to a shift, the sequence is symmetric. In Table I, we 
give the non-repeating strings of 15 admissible periodic sequences, the initial values (s,v) and the symmetry of their 
corresponding orbits. 

Orbits which pass through the origin play a special role in the semiclassical theory of closed orbits. Let us analyze 
their symmetry in some detail. Right at the origin the momentum v is ill-defined. However, for an orbit passing 
through the origin we may consider its left and right limits since both of them are well-defined. It can be seen 
that around the origin we have codes BDq and BD\ for the two limits, where B stands for Bo or B 2 . In the 
RD Dq should have been identified with D 2 since the dynamics is on the annulus. Here by Dq and D 2 we mean 
s = 0+ and 2_, respectively. It can be further verified that the first Dq(D 2 ) and second D 2 (D$) are related as 
(s, v) <-> (2 — s, v), which is just the transformation T. This means that, after identifying Dq with D 2 , the sequence of 
an orbit passing through the origin must be T-invariant. For example, the orbit shown in Fig. 9(c) has the sequences 
(B 2 D' 2 R B D 2 R 2 )°° and (BqD2R 2 B 2 D 2 Rq)°° as its right and left limits. Written more precisely, the two limiting 
sequences are (R 2 L 2 RiRqRoLqRiR 2 ) oc and (R 2 RiLqRqRqRiL 2 R 2 ) co , respectively. The T transformation of the 
former is (R R 1 L 2 R 2 R 2 RiL Ro) oa , which turns out to be just the 4-th shift of the latter. (A further observation 
discovers that the codes BDD ■ ■ ■ can be equally well written as BDB ■ ■ ■ since the third point is on the common 
segment of »D and »B, where D, similar to B, stands for D or D 2 .) 

Many short closed orbits have been found in Ref . pi . Non-repeating strings or their doubled strings of the periodic 
sequences corresponding to these orbits are listed in Table II, where sequences and initial values v are for the right 
limit orbits starting at s = + . Sequences for left limit orbits may be obtained by exchanging Dq with D 2 . While 
the labeling of Ref. || for orbits are given in the table, we arrange them according to the order of sequences and 
symmetry, and shall refer this numbering to orbits. More precisely speaking, the two sequences which relate to each 
other by a are always listed successively, and ordered according to one of the two without referring to the order of 
the other. The orbit shown in Fig. 9(c) is orbit 27. 

From the table it is clearly seen that symbolic coding provides a way to classify orbits. For example, according 
to the coding orbits 1, 3, 5, 15, 25 and 2, 4, 10, 20 form a big family of the type BD 2 R™BD 2 R™. When n is odd the 
non-repeating string of an orbit reduces to half. We have doubled it in the table to make it look more like that for an 
even n. There are two more orbits (29 and 30) for which doubled strings are given. We can see many pairs of orbits 
which should be created in a bifurcation related to »C. A pair of orbits 29 and 30 are given in Fig. 9(d) to show that 
they have the same topology. Other examples are orbits 6 and 8, orbits 11 and 13, and orbits 28 and 30. In the table 
8 orbits are obtained by the transformation a, which is indicated in the last column of the table. The initial values 
of v for two such orbits related by a are different only in their sign, so only one v is given. The sequences listed for 
these orbits then correspond to a left limit. 

VI. CONCLUSION 

In the above we have given a natural way to construct symbolic dynamics for the the diamagnetic Kepler problem. 
Our choice of the Poincare section avoids any ambiguity in identifying bounces. A lift of the phase space helps us to 
understand the dynamics. Tangencies between stable and unstable manifold foliations play an important role in the 
construction of symbolic dynamics. Based on the symbolic coding we have analyzed the rotational and reflectional 
symmetry, and the time reversal symmetry as well. Symbolic coding provides us a convenient way to classify orbits. 
For a billiard limit we have established the relation between our coding and the coding based on bounces. 
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Table I. Some orbits which do not pass through the origin. 



No. 


s 


V 


Sequence 


Period 


Symmetry 


Bound Codes 


1 


0.502590559992 


0.000000000377 


Ro Ri R\ Ri Ro 


6 


(T, T, T CT 




2 


0.525466094731 


0.083479979410 


RqLoRiRo 


6 




2 2 102 


3 


0.597726644479 


0.000000000105 


Ro 


4 


T CT 


2 4 


4 


0.627764854085 


0.000000000692 


p2 r>2 n»2 


6 


T CT 


2 3 02 2 


5 


0.634615582851 


0.343668565696 


R1L2R1R0 


4 


T CT 


212 


6 


0.634764050232 


0.000000004571 


p2 r>4 n»2 


8 


T CT 


2 3 0202 2 


7 


0.637932489171 


0.466976827386 


R1R2R1R0 


4 


O", T, T CT 


20 2 2 


8 


0.644118560693 


0.000000282982 


p2 r>12 r>2 


16 


T CT 


2 3 (02) 6 2 


9 


0.645908548293 


0.860542468032 


Ri i?2 L2R1L0 Ro Lo 


16 




1(02) 3 2 1201 


10 


0.653504728307 


0.022526445642 


R0L0R1 7?o 


26 




21(02) n 2 


11 


0.653578168108 


0.014059036967 


RoLoRj 6 Rq 


30 




21(02) 13 2 


12 


0.671107671805 


0.000000002635 


RoLoRiLqRo 


8 




210201 


13 


0.685136083296 


0.257524880819 


Ri R2 Ri Ro LqRi Rq 


16 




20 2 2 2 10(20) 3 2 
20 2 2 2 02 2 1 


14 


0.714849526728 


0.061364237283 


R1R2R1R0R1R0R1L2R1 


10 




15 


0.744712436526 


0.064776613634 


R\Lo Ro Ri R\ Ri 


8 




20120 3 
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Table II. Sequences of orbits passing through the origin. By Oi we mean the i-th orbit according to the numbering in the first column. 



No. 


label 121 


V 


Sequence 


Period 


Symmetry 


Bound Codes 


1 


30 


0.974971232297 


D0R1B2D0R1B0D0 


22 




21(02) 4 01(20) 4 
21(02) 4 10(20) 3 


2 


24 


0.972898880448 


DnRl BqDoRJ B D 


20 




3 


19 


0.970340776638 


d r\b 2 dIrIb d 


18 




21(02) 3 01(20) 3 

21(02) 3 10(20) 2 

21(02) 2 01(20) 2 

21(02) 2 2 12 2 020 

2 2 10 2 (20) 2 1(20) 2 

21(02) 2 1 3 020 


4 


13 


0.967086416031 


D RIB DIr\B D 


16 




5 


9 


0.962773942237 


DoRiB 2 DoRtB Do 


14 




6 


29 


0.959653990845 


Dq R\ R2 B2 D2 Ro R\ Bq B)() 


16 




7 


14 




D 2 RfR B Dl R 2 R\B 2 D 2 


16 


a Oq 


8 


28 


0.958697661319 


D R\L 2 B 2 DILqR\ Bo D 


16 




9 


15 




d 2 r{l b dIl 2 r\b 2 d 2 


16 


a Os 


201(20) 2 1 3 20 
21(02) 2 1020 
21(02) 2 210 2 20 
12(20) 2 1(20) 2 


10 


8 


0.956722093284 


D RfB D 2 RfB D 


12 




11 


27 


0.952489366877 


Dq R\ Rq Bq Dq R2 R\ Bq JDq 


14 




12 


23 




D 2 RiR 2 B 2 D 2 R RlB 2 D 2 


14 


a On 


13 


20 


0.950042745974 


DoR'iL B D 2 L 2 RlB Do 


14 




210201 3 20 


14 


22 




D 2 R\L 2 B 2 DlLoR\B 2 D 2 


14 


<r O13 


1 2 02012021 


15 


5 


0.947457897321 


d r\b 2 dIr\b d 


10 




21020120 


16 


17 


0.940875954726 


JDq Ri R2 B2 D2 Ro Ri Bq JDq 


12 




21020 2 12 2 


17 


18 




d 2 rIr b dIr 2 rIb 2 d 2 


12 


Oi6 


2 2 10 2 20120 


18 


12 


0.935006464896 


DoRlL 2 B 2 DlLoRlBoDo 


12 




21021 3 


19 


16 




D2R\ LqBq B)\ L*2R\ B2D2 


12 


<7 O Ol8 


201 3 201 


20 


3 


0.931024837844 


DoRiB D 2 R 1 B D 


8 




210210 


21 


11 


0.918172430278 


D Ri Ro B DqR 2 Ri B D 


10 




2102 2 10 2 


22 


10 




D 2 RxR 2 B 2 D\RoRxB 2 D 2 


10 


<j O21 


12 2 0120 2 


23 


7 


0.899732522729 


DoR\LoBoDlL 2 R x BoDo 


10 




2101 3 


24 


6 




D 2 R\ L 2 B 2 Dq L0R1 B 2 D 2 


10 


a O23 


1 2 0121 


25 


2 


0.891672484472 


DoB 2 D BoDo 


6 




2101 


26 


21 


0.864842403049 


DoR 2 2 B 2 D 2 oRlBoDo 


10 




210 3 12 2 


27 


4 


0.849666377360 


DoR 2 B 2 D 2 2 RoB Do 


8 




210 2 12 


28 


26 


0.804140836919 


DoR 2 RlL BoDf ) L 2 RlR BoDo 


20 




1 3 (20) 2 2 2 10 2 (20) 2 
210(02) 3 210(02) 3 
1 3 (20) 2 1 3 (20) 2 
l 2 


29 


25 


0.803121053916 


Do R 2 R'l Ro Bo D 2 R 2 R'l Ro Bo Do 


20 




30 


31 


0.790624022578 


Do L 2 R'l Lo Bo D 2 L 2 R'l Lo Bo Do 


20 




31 


1 


0.707106780832 


DoD'iDo 


4 
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FIG. 1. Boundary of the transformed potential for the diamagnetic Kepler problem at e = 0, and a typical orbit. 



FIG. 2. The image (a) and preimage (b) of the fundamental domain. 



FIG. 3. The partition of the reduced domain according to preimage (a) and to present (b). 

FIG. 4. (a) Symbolic plane of the reduced domain. Approximately 68000 real orbit points are shown, (b) Primary pruning 
front. Points on •D and »D 2 are marked with diamonds, and points on •Co and »C 2 with crosses. 

FIG. 5. Symbolic plane of the minimal plane. Approximately 34000 real orbit points are shown together with the primary 
pruning front. Points on »Dq are marked with diamonds, and points on »Co with crosses. 



FIG. 6. Admissible and forbidden orbits are determined by five points on two partition lines. 

FIG. 7. The touching four-disk billiard (a) and its phase space (b). Here the billiard boundary is taken as a Poincare section. 
The Birkhoff coordinates are denoted by (r,u). 

FIG. 8. The grazing orbit and center-passing orbit (a) determine the correspondence between (s, v) space and the 
bounce-based (r, u) space (b). 

FIG. 9. Examples of orbits: (a) u-symmetric (RiR^RiRq) 00 (solid) and asymmetric {R\Rq)° c (dashed); (b) an orbit 
(RlLoRoRiR 2 )°° (solid) and its TV-transformation image (RlR 1 R 2 L 2 Rl)° c (dashed); (c) an orbit passing through the ori- 
gin (orbit 27 in Table II); (d) orbits DoR 2 RiL B D^L 2 RIRoB D (solid) and D R 2 RIR B 2 DIR 2 RIRoB D (dashed) with 
the same topology. 
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